Ponderomotive acceleration of hot electrons in tenuous plasmas 



V. I. Geyko and G. M. Fraiman 
Institute of Applied Physics RAS, 46 Ulyanov St., Nizhny Novgorod, 603950, Russia 

I. Y. Dodin and N. J. Fisch 
Department of Astrophysical Sciences, Princeton University, Princeton, New Jersey 08544} USA 

(Dated: February 9, 2009) 

The oscillation-center Hamiltonian is derived for a relativistic electron injected with an arbitrary 
momentum in a linearly polarized laser pulse propagating in tenuous plasma, assuming that the 
pulse length is smaller than the plasma wavelength. For hot electrons generated at collisions with 
ions under intense laser drive, multiple regimes of ponderomotive acceleration are identified and the 
laser dispersion is shown to affect the process at plasma densities down to 10 17 cm -3 . Assuming 
a/y g <C f , which prevents net acceleration of the cold plasma, it is also shown that the normalized 
energy 7 of hot electrons accelerated from the initial energy 70 < V does not exceed V ~ ay g , where 
a is the normalized laser field, and 7 9 is the group velocity Lorentz factor. Yet 7 ~ F is attained 
within a wide range of initial conditions; hence a cutoff in the hot electron distribution is predicted. 
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I. INTRODUCTION 

The recent advances in the laser technology have 
yielded techniques for generating electromagnetic radi- 
ation with intensities as high as 10 22 W/cm 2 [1 . Experi- 
ments show that interaction of ultrapowerful pulses with 
underdense plasmas produce hot electrons with energies 
up to hundreds of Me V [H [5] . As argued in Ref . [I] , the 
effect might be due to ponderomotive acceleration of elec- 
trons, following large-angle collisions with ions in strong 
electromagnetic field. Assuming that the laser disper- 
sion is negligible due to the plasma density being small, 
the model explains the observed power-law spectra and 
predicts that the particle maximum energy scales as the 
third power of the field amplitude. This estimate is also 
in approximate agreement with the available experimen- 
tal data however, the latter is insufficient to conclude 
whether the model is, in fact, quantitatively accurate. 
On the other hand, already small yet nonvanishing den- 
sities of the plasma can undermine the assumption of 
negligible dispersion and therefore modify the acceler- 
ation mechanism: the electron velocity can then exceed 
the group velocity of a laser pulse, so the particles can be 
reflected, or "snow-plowed" by the field envelope. Thus, 
to understand the production of hot electrons in previous 
and future experiments, the effect of the laser dispersion 
on ponderomotive acceleration must be explored. 

Previously, the snow-plow acceleration was studied 
in specific regimes when the electron motion becomes 
exactly integrable. Particularly, Refs. [5J |7] as- 
sume equal group and phase velocities of the laser, and 
Refs. [8, 9, 10, llj suppose circular polarization and cold 
electrons (i.e., having zero transverse momentum), also 
adopted in Refs. [HI [T3J [TH Q~5] for an oscillation-center 
model. However a general treatment of the relativistic 
ponderomotive force in plasma has not been formulated, 
and the effect of the laser dispersion on the ponderomo- 
tive acceleration of hot particles has not been understood. 



The focus of this paper is then twofold. First, we de- 
rive the oscillation-center (OC) Hamiltonian for a rel- 
ativistic electron injected with an arbitrary momentum 
in a linearly polarized laser pulse propagating in tenuous 
plasma, assuming that the pulse length Ly is smaller than 
the plasma wavelength X p . Second, we use this formalism 
to describe the ponderomotive acceleration of hot elec- 
trons generated at collisions with ions under intense laser 
drive. Specifically, we identify multiple regimes of this ac- 
celeration and show that the laser dispersion affects the 
process at plasma densities down ton~ 10 17 cm -3 . As- 
suming a/jg <C 1, which prevents net acceleration of the 
cold plasma, we also show that the normalized energy 7 
of hot electrons accelerated from the initial energy 70 < T 
does not exceed T ~ a"f g , where a is the normalized laser 
field, and 7 9 is the group velocity Lorentz factor. Simul- 
taneously, 7 ~ r is attained in a wide range of initial 
conditions, with the angular spread of the accelerated 
electrons \ ~ Ig 1 - 

Hence we conclude that the distribution of hot elec- 
trons produced at large-angle collisions with ions at 
iy -C A p and a/~f g -C 1 must have a cutoff at the energy 
7 ~ aj g . This refines the result from Ref. [4], showing 
how even weak laser dispersion can affect the accelera- 
tion gain. However, further experiments are yet needed 
to validate the updated scaling, because no relevant data 
has been reported for the regime considered here. 

The paper is organized as follows. In Sec. [H] we in- 
troduce our basic equations. In Sec. |III| we derive the 
OC Hamiltonian for a particle interacting with a laser 



pulse in tenuous plasma. In Sec. IV we identify the ma- 
jor regimes of ponderomotive acceleration in plasma and 
find the general expression for the particle energy gain. 
In Sec. [V] we discuss what we call the plateau regime, 
where 7 ~ T is attained within a wide range of initial 
conditions. In Sec. I VII we summarize our main results. 
Supplementary calculations are given in Appendix. 
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II. BASIC EQUATIONS 

Suppose a plane laser wave propagating in plasma with 
the group velocity v g and the phase velocity v p along the 
x axis, so the vector potential reads A = y°A, 



A = A 



L 



v t 

'' ^ cos(k[x — v p t]) . 



(1) 



Here y° is a unit vector along the y axis, Ly is the spatial 
scale of the envelope A, and k is the wavenumber such 
that e = (/cZ/||) _1 <C 1. Consider a particle with mass m 
and charge e interacting with this wave, assuming that 
the electrostatic potential is negligible (Sec. VBl. Then 
the particle Hamiltonian is [16 



H = c\l m 2 c 2 



■n 



P± -A 

c 



(2) 



where p x is the x component of the particle kinetic mo- 
mentum, and P± is the conserved transverse canonical 
momentum. 

In the extended phase space, where (t,—H) serves as 
another canonical pair and the independent variable is 
the proper time r, the equivalent Hamiltonian reads [17 



8 = 



1 

2m 



Pi- 



ll 1 



(3) 



and, numerically, S) 
variables 



0. Introduce the dimensionless 



x — kx . 
i = ckt, 
f = kcr, 

f3g = Vg/C, 



P = Px/mc, 
7 = H/mc 2 , 
h = Sj/mc 2 , 
(3 P = Vp/c, 



(4a) 
(4b) 
(4c) 
(4d) 



and V±_ = Pj_ Jmc — const. Hence we rewrite Eq. ^ as 
1 



1 + p* + (9 y - A) 2 + 3* 



7 



(5) 



assuming r is the new time, and the normalized laser 
field A = eA/mc 2 reads 



A = a[e[x — 0gt]) cos(x — (3 p t). 



(0) 



III. OSCILLATION-CENTER HAMILTONIAN 

A. Extended Hamiltonian 

Like in Refs. [6l EJ [8j El [TUJ H2J [13l [15], we assume 
the linear plasma dispersion, which holds for arbitrarily 
large a at L n <C X p [T5J [HE [2U]. [Also, the nonlinear 
instabilities will be neglected as they occur on time scales 
exceeding the acc elerat ion time, which is less than the 
wave period (Sec. IV B I.] Then, 



f3 g = \/l - a, j3 p = l/Vl - a, 



(7) 



where a = n/n c , and n c is the critical density 121) . 
Perform a canonical transformation on Eq. ^ |22j : 



{x lP ;t,-i)^{e,v- t,w), 

governed by the generating function 

F=(x- f3 p t)V + (x- /3 g t)W. 

Then 

P = V + W, l=PpV + /3gW, 

so the new Hamiltonian reads 



(8) 
(9) 
(10) 




aV 2 
1 - a 



+ aw 2 + y 2 + 



+ [?„ - a(&/>) cos 9] 2 \, (11) 



and the new variables are given by 

7 - PgP 



1p = X-(3gt, W 



7 ~ P P P 



(12) 
(13) 



Unlike at f3 p = f3 g [6j |7], e.g., for vacuum (Ap- 
pendix [A]) , or the exactly integrable case of circular po- 
larization with zero Vy [5J GO QUI EI] > there are two inde- 
pendent coordinates 6 and ip entering Ti. here; hence we 
proceed as follows. Introduce the normalized momenta 

V = aT, W = aW, (14) 

which remain finite at a — *■ 0: hence the Hamiltonian 



n 



i 



v 2 



a 



1 - a 



+ W 2 + aT;+ 



+ a[J> y - a(e^) cos 0] 2 \. (15) 

Following the general perturbation theory [231 [2"4l 1231 [2"6"] . 
we now seek to map out the quiver dynamics. To do this, 
consider a canonical transformation 

(6,P; ^,W)^(e,n;^,W) (16) 

governed by the generating function 

f = en + ipw + s(9,n-,ii>,w). (17) 

Choose S such that fl and W are OC canonical momenta, 
i.e., the new Hamiltonian H(0, II; 'J, W) does not contain 
fast oscillations. Then 



1 



1 - a 



2nd e S+{d Sf 



0, (18) 
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the tilde standing for the quiver part, and 

f(i/>,6) = ia 2 («/0cos2<9 - 23^ a(eV>) cos 0. (19) 

At e <C 1, the terms containing d^S are negligible; 
thus, from Eq. (18), S is nearly independent of W, and 



* = tp + d w S w ip. 



(20) 



Hence Eq. ( 18 I rewrites as 



1 

I- a 



2nd e S+(dgS) 2 ~C 2 +a/tyi,0) = O, (21) 



where C 2 = ((dgS) 2 '), and the angular brackets denote 
averaging over 6*. Solving Eq. (21 1 yields 



S 



en+ dn 2 + C 2 + a(l-a)f{ip,9) d6, (22) 



where we chose the root which corresponds to V > 0, 

P = n + d e S. (23) 

Require that S does not contain a zeroth-order harmonic 
in 0; hence, due to Eqs. p0|, |22l, (p3l, C is found from 



2tt 



n 2 + C 2 + a(l-a)f(m,0) dQ = 2nn. (24) 



Sec. 


IIIC 


0.) T 


icn 



H : 



a{\ + y\ + 2$) + W 2 - 



*^(i-* 2 ), 



5 2 = 



2C 2 



(25) 

(26) 
(27) 



a(l — a)a 2 

Hence we integrate the motion in the variables (0, /7): 

n = const, 6 = r7r/(l - a) + const, (28) 
and the remaining canonical equations read 

W = -<9*H, V = W = W -d^S (29) 

B. Effective mass M* 

One can also revert to the space and time coordinates, 
which is done as follows. Apply the variable change 



n = an, W = aW, H = aH, 



(30) 



where H(0,/7; ^f,W) is the new Hamiltonian. Perform 
a canonical transformation 



(6,/7; *,W) -» (X,T X ; T, —IK) 



(31) 



governed by the generating function 



IK. 



Pp Pg Pp Pg 

Then = {p), IK = (7), and 



Pg 



Pp Pg 



(32) 



(33) 



Now return from the extended phase space to the phys- 
ical phase space, so that T becomes the independent vari- 
able. Hence the new Hamiltonian 



IK = y/M* + V 2 (34) 
is equivalent to that of a particle with an effective mass 

M, = VI + 2$, (35) 
where <I> = $(a, a, CP V , II) , 

n 



n = 



3~C Pg^x 



(36) 



VT=a 

is a constant determined by the initial conditions, and 
IP 2 = IP 2 , + IPj_ is the OC total momentum squared. Thus 
the average force on a particle due to the laser field, or 
the so-called ponderomotive force, reads 



(7) _1 VM„ 



(37) 



in the nonrclativistic case yielding F sa — V<&, where <!> is 
called the ponderomotive potential [271 12"51 [2^j |5TJ1 |5T] . 
From Eq. (26), it flows that the plasmon inertia de- 



creases the electron effective mass and <F As the treat- 
ment is expanded to arbitrary dispersion (other polariza- 
tions are allowed, too), it can also be shown that M* < M 
at P p > 1, and M* > M at (3 P < 1 in the general case. 
However, the sign of the square root in Eq. (22 1 (and fur- 
ther) must be chosen appropriately, accounting for the 
fact that V [Eq. ( 12 1] might then become negative. 



C. Explicit approximation for M* 

To find the Hamiltonian H and the effective mass ex- 



plicitly, solve for S using Eq. (24 1, which rewrites as 

r \/^ 2 + 9„)/2 + 1(0, q y ) S = 2nC\ (38) 
Jo 

where 

C = a^/n, q = y/a, f = .f/a 2 . (39) 
At £ <C 1, £q y <C 1, this yields an approximate solution 

(40) 



5 « Uql + 1/16. 
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1: <5 [Eq. (27)] vs. £ = aVa/II for different ^ = CP H /a 



FIG 

solid - numerical solution of Eq. (241 for £ < £*(q y ); dotted 
analytical approximation (40 1. Intersections with S 2 (£, r ,q±) 
[Eq. ( |44[ )] (dashed) yield the reflection points £ r (dots), and 
physically realizable are £ < £ r . Here CP Z = and a — > 0; 
hence £ r ~ £* for = 0, 0.5, 1, but £ r < for g H = 2. 
The exact £ r are close to those flowing from the analytical 
approximation, except at q y — 0.5 here, in which case the 
dotted and the dashed lines do not intersect. 



Then S <C 1, so the effective mass reads 

/ a 2 S 2 \ 
M. = M 1 



4M 2 



(41) 



where M = yl + a 2 /2 is the effective mass in vacuum 
[HI EU [32j [331 El [35j [36] . Particularly, for cold particles 
with y — (i.e., II w M), and a — a a /M (assuming 
relativistic modification of the critical density, with cvo = 
const), one gets = M[l — aoa 4 / (6AM 5 )], in agreement 
with Ref. [37]. 

Eq. pO] ) can also be extrapolated as follows. Eq. ((38|) 



must hold for any initial conditions; however, at large £, 
its right-hand side goes to zero, and on the left-hand side 
£ -2 becomes negligible in comparison with |/| > 1. On 
the other hand, the square root in Eq. ( 38 1 is supposed to 



remain positive and non vanishing due to the oscillating /. 
Thus there is no solution for 6 at £ > 1, meaning that 
there exists £*(q y ) such that any realizable £ satisfies 



(42) 



(see also Sec. HID). Yet the exact numerical solution of 
Eq. ( 38 1 for <5(£, q y ) and its domain is close to Eq. (|40l) 



for any £ from the interval [0,^*{q y )}, as seen in Fig. 



Therefore Eq. (40 1 roughly holds for any realizable £, anc 
Eq. (351 can be used to, at least, estimate M* explicitly. 



D. Reflection point 

Since £ < £*, a particle cannot enter a field with 
a > a* = ^Il/y^a; thus, if the maximum field exceeds 
particle is reflected. On the other hand, not all 



£ satisfying Eq. (42 1 can be physically realized; thus a 



particle may bounce off even weaker field. 

Specifically, the reflection condition is found from 



a 



a 2 (e*) 



(1 - S 2 ) 



iU 



(43) 



which is obtained using Eq. (251, together with H = 0. 
Suppose a <C a 2 ; then, at 'F = 0, being the condition 
of particle stopping in the frame traveling with the laser 
envelope, Eq. (43 1 yields 



S 2 (Zr,q±) = 1 + 2 (<?!-£" 



(44) 



for the reflection point £ r . Unlike £*(q y ), the value of £ r 
is then determined by both q y and q z ; hence £ r < £*, 
except at q z — and Iq^I < 1, for which case one can 
show £ r — > for a — > (Fig. [I) . 

With Eq. (40 1 used as an estimate for 5 1381 . one can 



further show that, in agreement with Refs. [6"ll?ll51l9"lllO|. 

£ r ~min{l, -7I 1 }, (45) 



assuming the inequality ( 42 ) . Thus reflection is impos 



sible at £ <C £ r and possible at £ ~ £ r , whereas larger £ 
cannot be realized. Therefore 



£ < £r < 1, 



(46) 



which also yields, from Eq. (40 1 and q±> q y , that 

<5 < 1- (47) 

IV. PONDEROMOTIVE ACCELERATION 

A. Basic equations 

The particle energy 7, as affected by the ponderomo- 
tive force (37 1, can now be calculated as follows. Use 



Eq. (10 1 together with Eos. (14 1 for V and W. Further, 
substitute V from Eq. ((231), with S found from Eq. (p2l 



and employ Eq. (29 1 for W, with 'F from Eq. (43 1; hence 



7 = -{\/n 2 + aa 2 5 2 /2 + a[(a 2 /2) cos 2(9 - 27 y a cos 9] ± ^(l - a) (n 2 - a[l + 9*_ + (a 2 /2)(l - <5 2 )]) }. (48) 

I 

Thus the energy retained outside the field is where the plus and the minus correspond to the particle 

overtaking the pulse and falling behind it, respectively. 

7 » = i{n ± ^(1 - a) (na-a[l + ??_])}, (49) 



5 



If no reflection occurs and the average-force approxi- 



mation [from which Eqs. (48), (49 1 are derived] holds on 
the time interval (— oo, +oo), then %<, matches the en- 
ergy before entering the field, due to the conservation of 



P±=\Pi 












I 


IV 


Too ~ n/a, 


(50) 


w « 2|pJ 


W f 


« 1 


w~l/(2 Pll ) 



w ~ apn/2 



unless II 2 ^> a(l + T 2 ^) and the particle is transmitted; 
otherwise Eq. ( 49 ) is Taylor-expanded as 



1 



7o 



n 2 + Pj 
2n 



(51) 



-1 



-1/2 



FIG. 2: Scalings for to [Eq. ( |53[ )] depending on the normalized 
initial momentum p (p z = 0, a < 1). The roman numbers 
tag distinct regimes. 



[cf . the exact solution ( A5 I for vacuum] . 
Hence can be found by substituting II from 

n 2 = w 2 ~aa 2 [f + S 2 /2]. 



(52) 



Here we employed Eqs. (23 1, ( |27| , (36 1, (39 1 and, using 
Eqs. ((To), Q, substituted P = Wl - a, with 



= 7o - PgPo 



(53) 



found from initial conditions (hence the index 0). If a 
particle is born inside the field (Sec. IV B |, the initial 6 
itself depends on n and must be found from Eqs. (27 1, 



(24 1 or, approximately, from Eq. (40 1; yet an estimate 
can be obtained as follows. From Eqs. (36 1, (39 1, one 
gets that n 2 /(aa 2 ) ~ £~ 2 > 1, the inequality being due 
to Eq. (|42]) . Together with Eq. (j47]) , this means that, for 
an estimate, the term proportional to 5 2 can be omitted 
in Eq. (52), and, since |/| ~ max{l, q y }, one finally gets 



n 2 ~ max{i 



aaTy}. 



(54) 



B. Regimes of hot electron acceleration 



Consider a hot electron produced inside a laser pulse, 



e.g., due to ionization or collision (Sec. IV C I, at some 6q 
and a of the order of the maximum amplitude a max . 
Hence, as the particle starts to oscillate, it attains 7 ~ 700 



already on a fraction of the oscillation period [Eq. (48 1], 
like described in Ref. [39]. To calculate the associated 
energy gain, suppose an initial momentum p = p /mc, 
for simplicity assuming p z = and a <C 1; thus 



= P± — acos#o ~ max{p^, a}, 



(55) 



whereas pn will denote the x-component of the particle 
kinetic momentum. Then one of the six regimes is real- 
ized, depending on how w [Eq. ( 53 1] is expanded (Fig. [2| , 
and more regimes appear due to Eqs. (54 1, (|55[) allowing 



multiple scalings for n and 3L, respectively. 

Below we limit our consideration to only a part of 
these regimes, because of the following. According to 
Eqs. (451, (46 1, a pulse with a maximum amplitude sat- 
> 1 will snow-plow cold electrons of the 
1 and ?_l = im 



isfying aa 2 
background plasma, which have n 



However, this would result in a significant electrostatic 
potential (ahead of the pulse) which is not included into 
the model; thus we assume 



aa,. 



«1, 



(56) 



so only few, hot electrons could be snow-plowed. Assum- 
ing also a»l, twelve distinct regimes persist (Fig. |3| , 
and those of primary interest are discussed below. 



Acceleration in vacuum 



Suppose that an electron is produced at rest, e.g., due 
to ionization [301 ST] ; then, from Eq. ( 55 ) , 



3 s j_ = ao = a cos ( 



(57) 



At a <^ 1, the pulse travels much faster than the particle; 
hence the weak dispersion due to plasma is inessential 
in this case. Then Eq. (52 1 yields II « w = 1, so f ~ 
a^a « 1 and II 2 > a(l 



-f 2 



j_ ) , both because of Eq. ( 56 1 . 
Therefore particle reflection from the pulse is impossible 
in this case (Sec. HID I, and Eq. (51 1 applies, yielding 



7oo = 1 + cio/2, 



(58) 



in agreement with Ref. |39j and regime I in Fig. [3] 

When a particle is born with positive 3> 1, stronger 
acceleration is predicted from Eq. ( |5T| ) due to reduced n. 
Indeed, suppose a small pitch angle \o ~ Pj_/Ph an d, 
again, neglect the plasma dispersion (a — > 0); then 



pi 



Mi 



< 1. 



Similarly, Eq. (51 1 holds, so one gets 



a 2 Pn 



Pi 



xIpV 



(59) 



(60) 



covering regimes III. 1 and IV. 1 in Fig. [3] Hence only a 
small fraction of electrons is accelerated efficiently, par- 



ticularly those with Xo S P\\ *C 1- 
imum energy now scales as 7 G 



than that flowing from Eq. (58 1 by the factor p K ^> 1 



a 2 p 



However, the max- 
which is bigger 
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w [Eq. ( 


53 


1] V y [Eq. ( 


55 


l] n [Eq. ( 


52 


1] 7so [Eq. ( 
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FIG. 3: Regimes of ponderomotive acceleration for electrons born inside laser field with initial momentum p. The dashed and 
dot-dashed lines separate different domains corresponding to the same energy 700 . The dot-dashed line is also a schematic of 
the curve (691, at which the particle velocity equals the pulse group velocity. Shaded is the "plateau" where the maximum 
energy ~ aa _1//2 is independent of p. The roman numbers are the same as for the corresponding domains in Fig. [2] 



Specifically, the described effect is anticipated at large- 
angle electron-ion collisions in tenuous plasmas [4]. Sup- 
pose a cold electron oscillating in a laser field with a 
quiver kinetic momentum ~ a and zero average ve- 
locity. (For the general case, see Fig. [3] and Sec. |V| ) 
Suppose further that this particle collides with an ion 
such that the momentum vector instantaneously rotates 
toward the pulse propagation direction, i.e., 



px ~ P~Xo- 



Then the maximum from Eq. ( 60 1 reads 



7oo ^ & 7 



(61) 



(62) 



the result being called the a 3 -effect |4j, and the angular 
spread of the accelerated electrons is 



D. Modification of the a 3 -effect in plasma 

Increasing the number of accelerated electrons requires 
higher plasma densities, and the a 3 -effect is modified in 
this case because of the laser dispersion; hence the en- 
ergy gain is calculated differently. Pa rticularly, for elec- 
trons with the initial conditions ( |6l| , one has w ~ a^ 1 
[Eq. (59 1; regime IV] and q y ~ 1 ; then Eq. (|54| yields 



'}. At a = aa^ <C 1 (regime IV. 1), 

,-1 



II 2 ~ maxjaa 2 , 

one obtains II ~ a -1 , so the reflection condition is not 
met, and the plasma effect is negligible. Suppose now 
that a > 1 (regime IV. 2). Then one gets 



n 



ova, 



(64) 



so it becomes possible to reflect electrons from the pulse, 
at least, for some 0q. (In vacuum, this effect is impossi- 
ble because particles could not travel faster than light.) 
Hence the final energy is estimated from Eq. ( 50 ) as 



X«3V7oo~a 2 <1. 



(63) 



aa 



-1/2 



(65) 



15*0.0 

FIG. 4: The final energy = p x /mc of electrons accelerated 
by a plane laser pulse in tenuous plasma vs. the normalized 
vector potential envelope a — eA/mc 2 at collision and cosine 
of the collision phase 6q\ a — 10 -3 , a max = 15. The elevation 
corresponds to the electrons being snow-plowed. 



and the angular spread of the accelerated electrons is 



X 



V±/loo ~ Vot < 1. 



Now rewrite Eq. ( 65 1 as 7, 



(66) 

~ 1 / 2 . Then a uni- 
form scaling is obtained, which covers both regimes IV. 1 
and IV. 2, accounting for how the et 3 -effect is modified 
with the plasma density: 



7o 



a x mm 



{!• 



r -V2 



(67) 



This agrees with the results of our numerical simulations. 
Specifically, at <r < 1 we observed the vacuum a 3 -effect, 
and electron reflection from a pulse was seen at 



> 4.4. 



(68) 



Hence a sharp dependence of 700 on whether particles are 
reflected or not [albeit the scaling holds for reflected and 
transmitted electrons equaly, as predicted from Eqs. (49 1, 
(50 1] and the abrupt elevation in Figs. [4] [5] both agreeing 
with Eqs. pBl ((67 1. 



FIG. 5: Same as in Fig. [4] for 60 = 0, with a max = 10. 
(a) a = 1CF 4 : numerical (solid) and analytical, 700 = a 3 , 
corresponding to the o 3 -effect (dashed); no particle reflection 
from the pulse, (b) a = 10~ 3 (numerical); the elevation cor- 
responds to the electrons being snow-plowed. 



which corresponds to the particle traveling at the pulse 
group velocity, with pj_ <C p\\ (dot-dashed in Fig. |3j. 
Hence the respective interactions are classified as follows. 

• In regimes III. 3 and IV. 2, a particle initially travels 
along x-axis slower than the pulse and is acceler- 
ated up to the energy (651. 



• In regime V.I, a particle initially travels along x- 
axis faster than the pulse. However, it gains addi- 
tional transverse momentum before it escapes from 
the field, resulting in the same energy gain (65 1. 



• In regime V.2, a particle is fast enough to run ahead 
of the pulse such that the energy 7 is not affected 
(700 ~ 7o) j a s opposed to vacuum where 700 ~ 70a 2 
would apply at arbitrarily large 70 (cf. IV. 1). 

• In all other regimes, the particle gains energy 
smaller than both 70 and that given by Eq. ( 65 I . 



Thus for an electron born inside a laser field one has 



7o 



< 



max{7o, T}, 



(70) 



V. PLATEAU REGIME 

A. Maximum energy gain 

Now consider a more realistic case when the electron 
is also preaccelerated by the pulse before the collision; 
hence we assume arbitrary initial conditions instead of 
Eq. (61 1. Similarly to Sec. IV D one can show that the 



acceleration is affected by plasma only in regimes III. 3, 
IV.2, V.I, and V.2 (Fig. [3). Those adjoin the curve 



Pn = 



pi 



(69) 



where T ~ ajg is the energy of a particle comoving with 
the pulse, with the transverse momentum p± ~a>l 
and the group velocity Lorentz factor 7 9 = a" 1 / 2 . 

Assuming 70 < T, the maximum (over #0) of the parti- 
cle final energy is then attained in the "plateau" formed 
by the domains III. 3, IV.2, V.I, where it is independent 
of the initial momentum p and so is the angular spread 
of the accelerated electrons: 



7o 



(71) 



Below we assess the feasibility of the plateau regime and 
suggest an estimate for the energy of hot electrons which 
can be produced in conceivable experiments. 
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B. Required parameters 

The one-dimensional (ID) model above neglects elec- 
tron escape from the accelerating field in the transverse 
direction. This is a valid approximation if 



Ar < kL x /p ± , 



(72) 



where At ~ A 1 !'/ 1 !' is the normalized proper time of the 
interaction, and A^f ~ 1 because the acceleration occurs 



on a single period (Sec. IV B). In the plateau regime, 
Eq. (43) yields 4> ~ a^fa; thus Eq. (72 1 rewrites as 



Lj_/X > (2tt a 



(73) 



where we took A for the laser wavelength, and p± ~ a. 
For narrower pulses, the energy gain would be somewhat 
lower than that predicted by Eq. (71 1, particularly for 
particles born at A -C a max , as also confirmed in our 
numerical simulations (Fig. [6]) . Nonetheless one can an- 
ticipate the ID scaling to hold for feasibly focused ul- 
traintense fields down to about a ~ 10~ 4 . Hence the 
laser dispersion should affect the electron acceleration at 
plasma densities down to about 10 17 cm~ 3 . 

Now let us estimate the influence of the previously 
neglected wake potential tp, which impedes the accel- 
eration because the associated electrostatic force is di- 
rected oppositely to the ponderomotive force [TO ] 142 ] 143] . 
The energy gain due to the electric field E v = — Wip 
is 7^ ~ eEpL/mc 2 , where L ~ k Af'Yao is the inter- 



action length, or kL ~ 7 2 . Assuming the wake spatial 
scale of about the plasma wavelength A p = (fc-y/a) -1 
and the density perturbation of the order of n, the Pois- 
son's equation gives a v = eE ip /(m(xu) ~ \[ol. Then 
7 V ~ a v kL ~ a' 1 / 2 , yielding 7^/700 ~ a -1 < 1, i.e., 
the wake is indeed negligible [44] . 

Hence Eq. (71 1 is a valid approximation for estimating 



the electron final energy For example, at laser intensity 
/ ~ 10 20 W/cm 2 and wavelength A ~ 1 /im, correspond- 
ing to a re [A/(l/xm)] [7/(1.37 x 10 18 W/cm 2 )] 1 / 2 - 10, 
and n ~ 10 17 cm -3 , corresponding to a ~ 10~ 4 , Eqs. (71) 



predict 7^ ~ 10 3 and x ~ 0.01. Therefore hot electrons 
can be accelerated to energies of a fraction of GeV and 
will be scattered within a small angle of 0.6°. 



400 



300 - 



^ 200 - 



100 





V f = 5000 






e = 500 ~~~~ 






/ ! = 225 






( = 50 







4 6 

a 
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FIG. 6: Same as in Fig. [3] for do = and a max = 10, but for 
a two-dimensional pulse a = a max exp(— x 2 /L 2 — y 2 /^) with 
L|l = 11A and different £ = L±/X, where A = 2n/k. 



Assuming a/j g <C 1 [Eq. ( 56 1 ] , which prevents net accel- 



eration of the cold plasma, we also show that the nor- 
malized energy 7 of electrons accelerated from the initial 



energy 70 < T does not exceed T 



where a is 



the normalized laser field, and 7 ff is the group velocity 
Lorcntz factor. Simultaneously, 7 ~ T is attained in a 
wide range of initial conditions, with the angular spread 
of the accelerated electrons x ~ lg 1 ■ Hence the distri- 
bution of hot electrons produced at large-angle collisions 
with ions at L K <C X p and a/^ g <C 1 will have a cutoff at 
7 ~ ajg. This refines the result from Ref. [4], showing 
how even weak laser dispersion can affect the accelera- 
tion gain. However, further experiments are yet needed 
to validate the updated scaling, because no relevant data 
has been reported for the regime considered here. 
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APPENDIX A: ENERGY GAIN IN VACUUM 



VI. CONCLUSIONS 

In this paper, we derive the oscillation-center Hamil- 
tonian for an electron injected with an arbitrary momen- 
tum in a linearly polarized laser pulse propagating in ten- 
uous plasma, assuming that the pulse length Ly is smaller 
than the plasma wavelength X p . We then use this formal- 
ism to describe the ponderomotive acceleration of hot 
electrons generated at collisions with ions under intense 
laser drive. Specifically, we identify multiple regimes of 
this acceleration and show that the laser dispersion af- 
fects the process at plasma densities down to 10 17 cm~ 3 . 



In the case of vacuum, a simplified solution is possible 
as follows. Perform a canonical transformation [53] 

(x,p;t,-j)^{e,p';rj,-w) (Al) 

using a generating function 

F = (x -t)p' -tw. (A2) 

Then 

6 = x — i , p = p, rj = t , w = 7 — p, (A3) 
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and the transformed extended Hamiltonian is given by 
ti = 1 - w 2 - 2wp' + + {Ty - a(eO) cos<9} 2 = 0. 
Then w is conserved, yielding an explicit solution for p 
1 



Hence the particle energy 7 = p + w is obtained, and 



outside the field one has [cf. Eq. (51)] 



1 + w 2 



2w 



(A5) 



2w 



1 - w 2 + ?l + {y y - a{ 



EV) COS 



(A4) 
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